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Abstract. We study collisions between pairs of bubbles nucleated in an ambient false vac¬ 
uum. For the first time, we include the effects of small initial (quantum) fluctuations around 
the instanton profiles describing the most likely initial bubble profile. Past studies of this 
problem neglect these fluctuations and work under the assumption that the collisions posess 
an exact SO(2,1) symmetry. We use three-dimensional lattice simulations to demonstrate 
that for double-well potentials, small initial perturbations to this symmetry can be amplified 
as the system evolves. Initially the amplification is well-described by linear perturbation 
theory around the SO(2,l) background, but the onset of strong nonlinearities amongst the 
fluctuations quickly leads to a drastic breaking of the original SO(2,1) symmetry and the 
production of oscillons in the collision region. We explore several single-field models, and 
we find it is hard to both realize inflation inside of a bubble and produce oscillons in a col¬ 
lision. Finally, we extend our results to a simple two-field model. The additional freedom 
allowed by the second field allows us to construct viable inflationary models that allow os- 
cillon production in collisions. The breaking of the SO(2,1) symmetry allows for a new class 
of observational signatures from bubble collisions that do not posess azimuthal symmetry, 
including the production of gravitational waves which cannot be supported by an SO(2,l) 
spacetime. 
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1 Introduction 

This is the third installment in our investigation of fluctuations around highly symmetric do¬ 
main wall collisions. We began with a linearized analysis in [1] and continued with numerical 
lattice simulations of colliding nearly planar walls in [2]. In this paper our focus is on the 
full nonlinear three-dimensional dynamics of collisions between pairs of bubbles nucleating 
in an ambient false vacuum. Such collisions are expected to occur in inflationary models 
where our observable universe is contained within one of these nucleated bubbles. Among 
current cosmological paradigms, this class of models are often discussed in the context of the 
multiverse, and more specifically the landscape [3, 4]. We explain more precisely our notion 
of the landscape and the embedding of inflationary bubble models below. Bubble collisions 
are also relevant to the study of first-order phase transitions occurring in both particle and 
condensed matter physics. In the cosmological context, the universe may have undergone a 
number of these transitions in its infancy. Of course, collisions between more than two bub¬ 
bles are required for the phase transition to complete, although individual collisions between 
pairs of bubbles also contribute. 
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Collisions between pairs of bubbles are most frequently considered in scenarios of open 
inflation following false vacuum decay. In this class of models, our local cosmological history 
proceeds as follows. A bubble nucleates within a surrounding false vacuum corresponding to 
the Big Bang “singularity” of our observable universe. The interior of the bubble undergoes 
A In a > 50 efolds of inflationary expansion to dilute the effects of spatial curvature and seed 
the density perturbations that will eventually collapse to form the structure we see in the 
universe today. Inflation inside the bubble ends and the universe (p)reheats, thus starting the 
standard thermal history of the universe. Meanwhile the scalar fields relax to a new deSitter 
false vacuum with p\ in order to explain the present accelerated expansion of 

the universe. Of course, since our universe arises from the nucleation of a bubble in this 
scenario, additional bubbles can nucleate in the ultra-large scale ambient false vacuum as 
well. This presents the interesting possibility of a collision between the bubble containing 
our observable universe and another bubble. 

The scenario described above is often invoked in the context of the landscape, and we 
now outline the prevailing picture for this embedding. Motivated by current ideas in high- 
energy theory, it is postulated that at the energy scales relevant for inflation there are many 
effective degrees of freedom beyond those of the Standard Model. For simplicity, consider the 
case where these degrees of freedom can be modelled as a large number of scalar quantum 
fields interacting through a potential and semiclassical general relativity.^ We also assume 
some mechanism exists to seed random initial conditions on the potential surface at different 
points in space. ^ The combination of a large number of dynamical scalar fields interacting 
through a potential with a mechanism to seed random initial conditions is what we will call 
the landscape paradigm. The fields then evolve through a combination of classical evolution 
and quantum effects such as stochastic inflation and quantum tunnelling. Some regions of the 
universe may become trapped in metastable minima with positive vacuum energy creating 
a large region of inflating false vacuum. Since the field can subsequently tunnel out of this 
false vacuum region through nucleation of a bubble, open inflation is naturally embedded 
into the global dynamics. 

In this paper, we restrict ourselves to collisions in either one- or two-field models. The 
single-field approximation to the landscape is common in the existing literature, and we 
demonstrate explicitly that the transition from single field models through the addition of 
even a single field can dramatically change the outcomes of collisions. Thus, we encourage 
some caution in drawing direct conclusions about bubble collisions on a very high-dimensional 
landscape. Furthermore, as mentioned above the full dynamics on a complicated potential 
will consist of much more than just bubble nucleations and in some regimes quantum in¬ 
terference effects may be important. We will not concern ourselves with these global issues 
and assume that in some region of the landscape the approximation of well-defined pairwise 
bubble collisions is adequate. 

Many studies have considered the phenomenology of open inflationary models arising 
from bubble nucleation [5-10]. Collisions between vacuum bubbles have also been considered, 
beginning with the work of Hawking, Moss and Stewart [11]. Early studies were motivated 

^This neglects various gauge and fermionic fields, as well as discrete degrees of freedom such as units of 
flux. However, it is sufficient to illustrate the setup. 

^As an example of such a seeding mechanism, some portion of the wavefunction for the (unknown) funda¬ 
mental degrees of freedom may describe QFT coupled to semiclassical general relativity. If this portion of the 
wavefunction decohered, we would naturally expect it to seed stochastic initial conditions for the emergent 
scalar fields. Dynamics on the potential surface (such as eternal inflation) may also serve to populate the 
landscape. 
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by the dynamics of early phase transitions within our horizon, with much of the focus on 
the production of gravitational waves [12-19]. Recent studies are more often motivated by 
bubble nucleation in false vacuum inflation [20-25]. These investigations have culminated 
in a general relativistic treatment of the bubble problem [26-29]. However, all of these 
studies assume the dynamics obeys an SO(2,l) symmetry, and no consideration is given to 
perturbations that break the symmetry. Even studies of fluctuations around the SO(2,l) 
collision restrict to perturbations that can still be evolved using the approximation of 1+1- 
dimensional field theory [30-32]. For the case of a single bubble the evolution of a more 
general set of fluctuations has been considered [33-35]. 

In an effort to observationally constrain such scenarios, several recent studies have 
proposed observational signatures from the collision between our “observation” bubble and 
an external “collision” bubble [36-40]. Several data searches have also been performed in 
conjunction with theoretical work, but thus far no signal has been found [41-47]. 

Since observing a collision with another bubble universe would revolutionize our under¬ 
standing of the cosmos, it is important to properly assess the possible outcomes of collisions. 
As well, given the fundamental nature of bubble nucleations to any first-order phase transi¬ 
tion, a thorough understanding of the collision dynamics is vital to a complete understanding 
of such transitions. Of particular relevance to us is the validity of the SO(2,l) symmetry as¬ 
sumption, which has never been explicitly checked. The results presented in this paper give 
a more complete understanding of bubble collisions in scalar field theories, without relying 
on the assumption of SO(2,l) symmetry to simplify the problem. We demonstrate that in 
a broad class of collisions, the dynamical amplification of quantum fluctuations leads to a 
complete breaking of the SO(2,l) symmetry. Since the SO(2,l) assumption severely con¬ 
strains the form of the final observable signals, it is possible that the searches performed 
thus far (which were based on studies assuming SO(2,l) symmetry) have been blind to other 
interesting possibilities. Although we focus on the dynamics of the collisions rather than the 
final observable signatures, our results suggest several interesting new observational avenues 
for testing open inflation. 

The remainder of the paper is organized as follows. In section 2 we present our models 
and outline a pseudospectral approach to solve for the initial bubble profile. Section 3 
considers the evolution of individual bubbles, both in the thin-wall and thick-wall regimes. 
A brief review of the collision dynamics under the assumption of SO(2,l) symmetry is given 
in section 4. We demonstrate that our lattice simulations reproduce the expected behaviour 
and test the sensitivity of the dynamics to the accuracy of the initial bubble profile. Our 
main results are then given in section 5 where we study collisions between bubbles in a 
variety of single-field models, including small initial fluctuations around the SO(2,l) profile. 
We extend these findings to a simple two-field model that permits inflation inside the bubble 
in section 6. Some brief comments on possible observable signatures and extensions to other 
scenarios are given in sections 7 and 8. Finally we conclude in section 9. 

2 Single Field Models and Description of Initial Conditions 

To provide a concrete setting we focus primarily on double-well potentials with broken Z 2 
symmetry.^ We make two choices for the symmetry breaking 



( 2 . 1 ) 


^Additional potentials are introduced later in the text as needed. 
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Figure 1. Plots of the symmetry breaking potentials for several choices of S. On the left we show the 
linear symmetry breaking (2.1) and on the the right the cubic symmetry breaking (2.2). The largest 
value of S in the linear potential is that for which the second minimum disappears. 


and 

Vcubic((/>) ^ ^ ~ ^o) + ^^^0 (^^3 ~ ^ 3 ) 

which we refer to as the linear and cubic symmetry breaking potentials respectively. In both 
cases S controls the difference between the false and true vacuum energies and Vq can be 
adjusted to give the desired true vacuum energy. Our investigations here are restricted to 
Minkowski space, so Vq does not influence the dynamics.^ In order for solutions describing 
vacuum tunnelling to exist in Minkowski spacetime, the false vacuum and true vacuum 
must be nondegenerate. For the linear potential (2.1), the second minimum disappears 
for |d| > 2/\/3, while for the cubic potential (2.2) this occurs when |d| > 1. In the linear 
potential, the locations of the two minima depend on d, while they are fixed at ±00 for the 
cubic potential. Unless explicitly stated, dimensionful parameters are measured in units of 
m = V^0o with the exception of the field measured in units of 0o Ihe energy density p 
measured in units of A0 q. 

2.1 Solution of the Instanton Equation 

Before solving for the dynamics of bubble collisions, we must first find the initial conditions 
describing the profiles of the bubbles at nucleation. In this subsection we present a new 
and extremely accurate numerical approach to determine the shape of the nucleated bubbles. 
Throughout this paper we assume the most likely bubble to nucleate within a surrounding 
false vacuum possesses an SO {A) symmetry in Euclidean time and is described using the 
bounce formalism of Coleman [48-50]. In Minkowski space, the profile satisfies the Euclidean 
signature equation 

^ + (2 3) 

drl ^ VE dvE 90 “ ^ ’ 

with the boundary conditions 

0(rE = 00 ) = 0false = 0) = 0. (2.4) 

^We also performed several runs with a homogeneous background Hubble and found no qualitative change 
to our results. 
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Here + (ix^ is the Euclidean radius and r = it is the Euclidean time. For a single 

scalar field in Minkowski space with potential satisfying some regularity conditions, it has 
been shown that the solution to (2.3) is the minimum action solution relevant to false vacuum 
decay [51]. 

Before presenting our technique, we briefly review the most common analytic approach 
(the thin-wall approximation) and numerical approach (the overshoot-undershoot method) 
used to solve (2.3). In the thin-wall limit, valid if the initial radius of the bubble is much 
greater than the width of the wall, we can obtain simple analytic estimates for the bounce 
profile. For the potentials (2.1) and (2.2), this approximation is valid when 5 <C 1. In the 
simplest form of the approximation, we drop the friction term and the terms proportional to 
S in the equation of motion. This leaves the equation for a domain wall in the degenerate 
double well ^ (0^ — (/)q) . Performing the usual quadrature gives an expression for the field 
profile 


= 00 f anh 


r - i?o 
V2 


(2.5) 


We then adjust the multiplier on the tanh and add a constant so that the field interpolates 
between the false and true vacuum. Finally, conservation of energy gives the initial bubble 
radius as 

i?o = — where a = 

Ap 


/ 


dr[dr(f)waii{i = 0 )]" 


(2.6) 


is the surface tension of the bubble wall and Ap = V{(j)false) — Vi^^true) is the difference in 
potential energy between the false and true vacuum. For the linear potential (2.1) we have 


^linear ^ ( 2 . 7 ) 

and for the cubic potential (2.2) we have 

v2 


It should be clear that with some modifications the thin-wall approximation can be applied 
to other potentials. In particular, if we replace 4)true with the location where the field tunnels 
out {(j^tunnel)^ wc Can apply the thin-wall approximation in situations where the field does 
not tunnel out near the true vacuum. Of course, we must either estimate <j)tunnel oi* ^Ise 
determine it by solving the bounce equation (2.3) in order to obtain the initial radius. If the 
latter approach is used to find (^tunneh then the thin-wall approximation itself will provide 
little utility other than aiding intuition. 

Although the simplicity of the final result makes the thin-wall approximation attractive, 
there are situations where it is either invalid or else insufficiently accurate. We must then 
turn to a numerical solution of (2.3). In the literature, the equation is usually solved via a 
shooting method. The ODE is recast into a root finding problem 0(r^ = (X)|0(r^ = 0) = 
0o) ~ ^ false — 0 for the initial field value 0o at te = 0. For the single-field case, the root 
finding is often done via bisection and the resulting algorithm is known as the overshoot- 
undershoot method. To extend this to the multifield case one needs to use a root-finder 
that works in multiple dimensions. An obvious choice is a Newton-type method, with the 
required derivatives obtained by solving for a collection of nearby trajectories and then doing 
a polynomial interpolation. Of course, the desired solution is a saddle point and care must be 
taken to evaluate these derivatives before the neighbouring trajectories diverge away from the 


- 5 - 




10° 


- 5 = 0.005, lin 

1 1 1 

- S = 0.001, lin 



- S = 0.01, cub 

10-4 


— (5 = 0.005, cub 

^ 10-8 


- 

10-12 


if™. 

10-16 

1 1 1 


) 20 40 60 

80 100 120 140 


Mode Number (i) 


Figure 2. Numerical solutions to (2.3) in the thin-wall limit. For each case the ratio of the wall width 
{rnw = \/2) to the initial bubble radius is w/Rq ^ 5. In the left panel we show the instanton profiles 
(zoomed into the region ve ^ Ro), while in the right panel we show the resulting spectral coefficients. 
For all cases, we took L = V^Rq with Rq = for the linear potential and Rq = 1.5\/2S~^ for the 

cubic potential. Meanwhile, the “stretching” parameter d was d = 0.022,0.0045 for 6 = 0.005,0.001 
in the linear potential, and d = 0.03, 0.014 for 6 = 0.01, 0.005 in the cubic potential. 


target solution. One possible way to extend the overshoot-undershoot approach to multiple 
fields can be found in [52]. 

Rather than adopt a shooting method, we instead use a pseudospectral approach and 
expand the function in even rational Chebyshev functions on the interval (—oo, oo) 


4^bounceij' e) ^ 


2i y 


y{x) — — tan ( d ^ tan ( tt 



(2.9) 


where are the Chebyshev polynomials. Since this is a global expansion method, it displays 
exponential convergence as the number of lattice sites is increased and machine-precision ac¬ 
curacy is easily obtained. Without the mapping y{x) (i.e. setting y{x) = x), this is simply an 

expansion in the rational Chebyshev functions TBn{y) = B^ ~ {ncot~^ (l)) 

on the doubly-infinite interval. L is an adjustable parameter that determines where the 
oscillations in the rational Chebyshev functions (or equivalently the collocation points) are 
clustered on the infinite interval. The even TB{x)^s form a complete set for smooth functions 
which asymptote to a constant at infinity and are symmetric about the origin (such as the 
bounce). Thus they enforce the = 0 boundary condition automatically. We also include 
an additional mapping y{x) which clusters (repels) points around if d < 1 (d > 1) 

without simultaneously introducing new singularities at the boundaries which would greatly 
reduce the convergence rate of the expansion. Via a judicious choice of L and d, we can 
resolve instantons to machine precision even in the extremely thin-wall case. We illustrate 
the efficiency of this expansion in figure 2, where we show instantons for four extremely thin- 
walled bubbles as well as the spectral coefficients q in (2.9). Numerically this is the most 
difficult case if no finesse is used in the solution. In all cases the round-off plateau for the 
coefficients is clearly present at i 100, indicating we have hit the limits of double-precision 
arithmetic. Through the use of a smart collocation grid, we manage to achieve this accuracy 
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Figure 3. In the top row we show instanton profiles in the linearly asymmetric double well (2.1) for 
a range of S values. The left panel shows the instanton profiles (zoomed into the region where the 
field is varying). Meanwhile the right panel shows (/)(rE = 0) which becomes the initial field value at 
the center of the bubble. For comparison we also include the location of the true vacuum minimum. 
The bottom row shows the same two plots for the cubically broken symmetry (2.2). For S > 0.25, 0.4 
the location where the field tunnels out to begins to deviate significantly from the false vacuum. 


using fewer modes than the ratio of the bubble radius to its width As an added bonus, 

our outer collocation point is located at = oo so there are no errors introduced by putting 
the system in a finite box. To demonstrate the utility of this approach across a range of 
potential deformations, figure 3 shows the instanton profiles in the linear (2.1) and cubic 
potential (2.2) for a range of S. 

Of course, the bounce solution only describes the most likely initial profile for the field. 
There are additional fluctuations which we include by taking the initial conditions for a single 
bubble to be 

0zmt(x, 0) = (t)bounce{rE = |x|) -h 50 0) = 50 (2.10) 

where (^bounce is a solution of (2.3) and 50 and 50 are realizations of random fields. This 
is generalized to multibubble initial conditions in the obvious way. The fluctuations 50 
encapsulate the effects of quantum fluctuations: deviations of the nucleated bubble wall 
from perfect SO(3,l) symmetry, bulk fluctuations in the ambient false vacuum, and bulk 
fluctuations in the bubble interior. Since we consider scenarios where the bubbles form via 
quantum nucleation, it is inconsistent to ignore the fluctuations since rare coherent excursions 
of these fluctuations are what allows nucleation to occur at all. 

A proper determination of the initial conditions requires a calculation of how the fluctu¬ 
ations in the original false vacuum are projected into the bubble spacetime by the nucleation 
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event (see e.g. [34, 35, 53]). This calculation is beyond the scope of this paper and would 
only serve to obscure the essence of our result. Instead we simply initialize the fluctuations 
as a realization of a homogeneous Gaussian random field with spectrum 


{\4k? 




{\4k? 




( 2 . 11 ) 


to mimic the fluctuations if the field were sitting at its false vacuum mimimum in Minkowski 
space. ^ The overall scale of the potential A enters into the initial fluctuation amplitude as 

^ (X V^. We measure the fields in units of 00 and time in units of ? so the 

vev 00 only appears in the equations of motion if we consider coupling to gravity where it 
determines the strength of the gravitational interaction via We use the same convolution 
based method as DEFROST to initialize the fluctuations [54, 55]. 

The remainder of the paper presents results obtained from a massively parallel scalar 
field lattice code written by one of the authors. Hamilton’s equations for the discretized 
system were evolved using a sixth-order symplectic Yoshida integrator [56-58] and (unless 
otherwise indicated) a second-order accurate and fourth-order isotropic stencil for the Lapla- 
cian [54, 59]. All production runs used lattices with 1024^ sites, although we did perform 
some numerical checks using 2048^ lattices. For all cases, the total energy of the system 
(when running Minkowski simulations) was conserved to the 10“^ level or better. 

To provide a test of our lattice simulations and facilitate comparison with previous 
studies, we also present the results of dimensionally reduced 1-h 1-dimensional simulations at 
several points in the subsequent sections. These simulations used a Fourier pseudospectral 
lattice discretization and a 10th order Gauss-Legendre time integrator. The basics of this 
approach are outlined in appendix B of [1]. 


3 Evolution of a Single Bubble 

First we consider the evolution of individual bubbles. Fluctuations in the angular dependence 
of the bubble radius for the case of individual bubbles in the thin-wall limit were studied and 
found to be stable in [33, 35]. Our results in this section demonstrate that this is also true 
for the bulk fluctuations around single bubbles. We evolve a thin-walled bubble with 6 = 0.1 
in the linear potential and a thick-walled bubble with 6 = 0.99 in the cubic potential. The 
corresponding instanton profiles can be found in the left panels of figure 3. For the thin- 
walled bubble, the field tunnels out very close to the true vacuum. As a result, we can think 
of this as a bubble with true vacuum interior separated from the false vacuum exterior by a 
domain wall of width m~^. For the single bubble this wall can be approximated as infinitely 
thin. The pressure differential between the interior and exterior of the bubble then causes the 
wall to accelerate and it follows a hyperbolic curve given by r^all{t) — ^/'^wall{i = 0)^ +as 
shown in figure 4. 

For the thick-walled bubble the evolution is somewhat different as seen in figure 5. The 
field now tunnels out far from the true vacuum. In the bubble interior it begins to oscillate 

^As a check of this assumption, we also ran simulations with several other choices of initial fluctuations. 
One set included homogeneous bulk fluctuations with different initial spectra than the Minkowski vacuum. 
Another set involved initializing 0init = 0bounce(r£; = |x| + 5r) with 5r — a 2d random field 

obtained by realizing a collection of a^rn’s. Again, we used several different spectra for our a^rn’s to verify 
that our results were not sensitive to a particular choice. In every case we tried, the outcome of a collision of 
two bubbles was qualitatively the same as the results we present in the remainder of the paper. 









mx_\_ 


mx 


Figure 4. Time evolution of a thin-walled bubble with 6 = 0.1 in the linear symmetry breaking 
potential. In the left panel we show the field value as a function of mt and position mx±_ along a slice 
taken through the center of the bubble. This figure is from a three-dimensional lattice simulation. 
For reference, the location of the wall in the thin-wall approximation is included as a dashed black 
line. In the right panel, we show field profiles at time mt = 0,17.24,34.48 on a slice through the 
center of the bubble. To demonstrate the accuracy of our numerics, we compare our lattice results 
to higher-resolution one-dimensional simulations. For comparison, we also include the evolution of a 
bubble whose initial profile is given by the thin-wall approximation rather than an exact numerical 
calculation. 


around the minimum at 0 = In the standard slicing of Minkowski space used in the 
code, these oscillations appear as outward propagating spherical waves. The leading edge 
of this spherical wave quickly develops into the bubble wall and propagates outward with a 
speed asymptotically approaching the speed of light. This is easiest to see by foliating the 
spacetime with hyperboloids centered on the nucleation site of the bubble. We label each 
of these hyperboloids by a new time coordinate t related to the time coordinate used in our 
simulations via t — tcoshx where x is the radial coordinate along the hyperboloids. The line 
element in the new slicing is ds^ — —dP + PdH^ with the unit three-hyperboloid. For 
the exact instanton initial condition, the evolved held is a function of i only and satishes 




(3.1) 


In this particular case the oscillating held in the interior of the bubble does not lead to a 
strong preheating instability because the potential seen by the held as it oscillates in the 
interior of the bubble is very nearly quadratic.^ To see this explicitly, dehne = 0 — 00 and 
reexpand the potential to obtain 


Uii^) = VicPo + ^) = (A + . 

The equation for linear perturbations around '0(t) is then 


(3.2) 


d{m^ffiy 


+ 




+ (1 + 5) + (3 + 5)0 + - ^ — 0 


2 


(3.3) 


®For different choices of potential or couplings to additional fields the interior of the bubble may experience 
strong instabilities (see for e.g. [60]). 
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Figure 5. Time evolution of a single thick-walled bubble with 6 = 0.99 in the cubic symmetry 
breaking potential. In the top panel, we plot the field as a function of mt and position mx±_ along a 
line through the center of the bubble. Regions where the field is near the false vacuum are blue, while 
regions where the field is near the true vacuum are red. In the bottom left panel, we plot the value of 
the field at the center of the bubble as a function of time. As expected the field undergoes damped 
harmonic oscillations. As a test of our numerical code, we include the result for two different choices 
of lattice spacing dx, and also for a one-dimensional simulation using a pseudospectral discretization 
and 10th order Gauss-Legendre time integrator. In the bottom right panel we show the field profile 
at times mt = 0,22.98 and 41.37, again for the two different lattice grid spacings and also the one¬ 
dimensional simulation. 

where = 0) = 2 A(/)q and k? is an eigenvalue of the Laplacian on the unit 

three-hyperboloid. To gain some intuition about these instabilities, let’s approximate the 
motion of ^ as it oscillates around the minimum hy ^ sin(me//t + do) with a < 1 

a numerical coefficient and Oq a phase. At the center of the bubble the time coordinate i 
coincides with the time parameter t used in our simulations. Therefore, we see from the 
middle panel of figure 5 that this parameterization of is indeed a good approximation. 
Ignoring the time dependence of all coefficients, the fluctuations obey an equation of the 
form 

dttf + + (3 -I- S)Bsm{P) + sin^(t)^ / = 0 . (3.4) 

For the continuum part of the spectrum, we have > 1 so A > (1 + d) and from figure 5 
we see that B <0.5 during the oscillations. In figure 6 we see that this puts us well into the 
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Figure 6. Floquet chart for / + (^4 + 4B sm{t) + 1.55^ sin^(t)) / = 0, corresponding to oscillations 
about the true vacuum minimum for (5 = 1. For oscillations in the interior of the thick-walled bubble, 
we have A>2 and B < 0.5. 

weak resonance regime. Since the oscillations damp with time, a given mode will trace a line 
in the instability chart so that even in the long-time limit there is no exponential growth and 
the decay is perturbative [61]. 

4 3D Bubble Collisions Without Bulk Fluctuations 

In this section we study collisions between pairs of nucleated vacuum bubbles working under 
the assumption of SO(2,l) symmetry. This allows us to compare our results with previous 
treatments of this problem (where the SO(2,l) symmetry is assumed to hold exactly and the 
dynamics are reduced to one spatial dimension). These simulations also provide a highly non¬ 
trivial test of our numerical approach. We consider both the case of the spectrally accurate 
numerical profiles for the bubbles, as well as initial profiles given by the thin-wall approxi¬ 
mation. Since the thin-wall profile is not an exact solution to the bounce equation (2.3), the 
latter choice is equivalent to a small breaking of the boost symmetry for the single bubbles. 
As we demonstrate below, even this minimal breaking of the symmetries of the bubble can 
have a dramatic effect on the full three-dimensional evolution. The nonlinear field dynamics 
results in a complete breakdown of the boost symmetries shortly after the collision. This 
breaking of the boost symmetries is our first new result requiring three-dimensional (more 
precisely higher than one-dimensional) simulations. 

4.1 Evolution of SO(2,l) Background and Linear Evolution of Fluctuations 

Before attacking the full problem, let’s review linear stability analysis around the SO(2,l) 
symmetric solution [1]. For the case of perturbations to a single bubble see [33-35]. To make 
the symmetry manifest, define new coordinates 

t = s cosh X 

X — X 

y — s sinh y cos 6 

2 ; = s sinhysin^ . (4.1) 
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Field configurations that depend only on s and x are SO(2,l) symmetric. Now consider 
fluctuations around such a symmetric background 


0 = 4>bgis, X) + 5(/>(s, X, X, (t>) ■ 


(4.2) 


In the limit of linear fluctuations, the background solution satisfies 

d‘^(t>bq , 2d(t>bq 


ds‘^ 


+ - 


s ds 


dx‘^ 


+ V'{c^tg)=0 


(4.3) 


and the fluctuations obey a one dimensional wave-equation with time- and space-dependent 
effective mass 


d^Ai 2 dAi d^Ai 


s ds 


dx‘^ 


+ \^^ + v" 


Af = 0. 


We expanded the perturbation in eigenmodes 


50 = A£{s^ x)C£^n{x)^^^^ with n G Z . 

i,n 

The Ci n labelled by I are the eigenfunctions and eigenvalues of 


d 


sinh(x) dx 


smhU)f) = + ^ 


c. 


(4.4) 


(4.5) 


(4.6) 


Past studies based on exact SO(2,l) symmetry restrict the treatment to a study of (4.3). A 
sample collision between two bubbles in the linear symmetry breaking potential with 5 = 0.1 
and SO(2,l) symmetry imposed is shown in figure 7. The bubble walls do not immediately 


£ 



mx 


o 



Figure 7. SO(2,1) symmetric collision between two bubbles in the linear symmetry breaking potential 
with S = 0.1. Red corresponds to regions where the field is near the true vacuum minimum and blue 
to regions where it is near the false vacuum minimum. The two bubble walls do not immediately 
annihilate each other in the collision, but instead they undergo multiple bounces while slowly radiating 
energy into the bulk. 

annihilate each other upon colliding, but instead bounce off of each other multiple times. Each 
bounce opens up a pocket at the collision site where the field is localized near the false vacuum 
minimum. The bouncing behaviour was first noted by Hawking, Moss and Stewart [11] and 
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is typical of thin-wall bubble collisions in double-well potentials. The bouncing behaviour 
can persist when gravitational interactions are included [26]. Due to the SO(2,l) symmetry, 
each pocket corresponds to an expanding torus in the full 3-dimensional evolution, where the 
field inside the torus is near the the false vacuum minimum. In [2] we showed that long tubes 
with the field near the false vacuum in the interior (formed from the breakup of colliding 
domain walls) are unstable to collapse. Hence, we expect that when the cylindrical symmetry 
is not perfect the torii will fracture into oscillons. In addition, small amplitude fluctuations 
localized near the axis connecting the bubble centers will see an oscillating background field 
as the two bubble walls bounce off each other. As a result, the fluctuations in these regions 
will be resonantly amplified, analogous to the situation with planar domain walls. We will see 
in section 5 that both of these behaviours manifest themselves in the fully 3-h 1-dimensional 
problem. 

4.2 3D Simulation of SO(2,l) Bubble Collsions 

Before proceeding to the case with bulk fluctuations around the pair of nucleated bubbles, we 
first present results using initial conditions that preserve the SO(2,l) symmetry up to errors 
induced from using a superposition of single instanton solutions. In the absence of numerical 
errors, the resulting time-evolution preserves the full SO(2,l) symmetry. Therefore, this is 
a highly nontrivial test of our numerics and also provides a nice visualization of the three- 
dimensional field profile. 

We consider our fiducial thin-walled instanton with 5 = 0.1 in the linear symmetry 
breaking potential (2.1). The single instanton profile can be seen in the top left panel 
of figure 3. From figure 8 we see that the lattice preserves both the rotational and boost 
symmetries quite well, especially considering that spatial discretization and discrete time 
steps explicitly break both symmetries. For times mt > 80 the surfaces of constant field 
are deformed slightly from being perfect hyperboloids near the edges of the collision, but 
they still exhibit the correct qualitative form. A late time three-dimensional view of the 
same collision as figure 7 can be found in figure 9, where the preservation of the rotational 
symmetry about the collision axis is clear. The expanding torii formed from the excursions 
of the field back to the false vacuum in the collision region are also clearly visible. The total 
energy of the system is conserved to the 10“^ level, and the deviation from perfect boost 
symmetry persists when we halve the time step or change the order of the time integrator. 
This indicates that the mild breaking of the boost symmetry is an artifact of the modified 
dispersion relationship of the finite-difference discretization rather than the time stepping. 

4.3 Sensitivity of Boost Symmetry to Initial Conditions 

We now consider the sensitivity of our results to initial perturbations by breaking the boost 
symmetries while still preserving the rotational symmetry about the collision axis. Specifi¬ 
cally, we perturb the instanton solution used to set initial conditions while assuming it is still 
only a function of r^. The resulting instanton profile no longer satisfies (2.3), and therefore 
the resulting evolution of a single bubble is not boost invariant. However, the three rotational 
symmetries are preserved. To see that the boost symmetries are broken, suppose that we 
choose a single bubble initial field configuration 0zmt,sym(x, 0) such that the time-evolved 
field is a function of |xp — alone (i.e. it is SO(3,l) invariant). Upon Wick rotation into Eu¬ 
clidean time, the SO(3,l) invariance translates into an SO(4) invariance so that the resulting 
Euclidean profile is a function of alone. As well, since the time evolved profile is obtained 
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Figure 8. The evolution of (j) along the line x = 0 = y where we have chosen our coordinates such that 
the two bubbles first collide at the origin and the collision occurs along the x-axis. For reference, we 
also include several lines of constant s, demonstrating that the boost symmetry is quite well preserved 
in the full three-dimensional simulation. 



Tlme=91.9239 Time=91.9239 


Figure 9. The post-collision field configuration for a two bubble collision in the linear symmetry 
breaking potential with 6 = 0.1 and exact initial bubble profiles. Left: Two slices of the field 
at mt = 92. The vertical slice is taken along the collision axis and cuts through the centers of 
the bubbles, while the horizontal slice is taken orthogonal to the collision axis at the collision site. 
Right: Contours of the energy density showing the expanding concentric torii whose interiors 
are near the false vacuum and whose walls interpolate between the false vacuum interior and the true 
vacuum exterior. An animation for the evolving field profile can be found at www.cita.utoronto. 
ca/~jbraden/Movies/linear_delO.l_exactprofile_field.avi. 


by solving the Klein-Gordon equation in Lorentzian signature, the Euclidean profile must sat¬ 
isfy the Euclidean version of the Klein-Gordon equation. However, for a function of alone 
the Euclidean Klein-Gordon equation is simply the equation for the bounce (2.3).^ Since our 

^This argument does not rely on the boundary conditions for the bounce solutions being met, only that it 
satisfies the correct equation. 
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perturbed initial condition does not satisfy (2.3), the resulting evolved bubble profile cannot 
respect SO(3,l) symmetry. By construction the rotational symmetry is preserved. Therefore, 
the use of an approximate initial field profile must destroy the boost symmetries. An example 
of exactly this type of approximate intitial condition is to use the thin-wall approximation 
rather than the exact instanton solution. Since the thin-wall approximation does not actually 
satisfy the instanton equation, this means it cannot respect the SO(3,l) symmetry. When we 
obtain two bubble initial conditions by superposing two of the approximate initial conditions, 
the evolved fields will preserve the rotation symmetry about the collision axis but break the 
two boost symmetries. 

Here we consider the effect of using the thin-wall approximation combined with an 
approximate determination of (pMse (/)true- This is a common approximation and is 
obtained if one writes the potential as Vsym + AH and solves for the profile of the wall by 
dropping the AH term (provided of course that the locations of the minima are perturbed 
by AH). Since we want the bubble profile to interpolate between the false and true vacua, 
we take the following initial condition for a bubble centered at xq 


(pinit = ^ tanh 


|x — Xo| — i?' 


init) 


W 


+ 


+ ' 


(4.7) 


with mi?init = and mw = y/2. The notation is meant to distinguish the ap¬ 

proximate locations of the false and true vacua from the exact values false/true- We further 
approximate the values of and to linear order in 5. We find 0^ = —l-h5/2-h 0{6^) 
and = 1 -h 5/2 -h (9(5^). At this level of approximation, we have 

^init ~ 00 tanh (li^^—-I- ^ . (4.8) 

A comparison of approximation (4.8), the thin-wall approximation with an exact determina¬ 
tion of (i>false/true the numerical result are shown in figure 10. Clearly, the approximate 

solutions provide a very accurate description of the system initially, although the perturba¬ 
tion to the initial bubble radius is visible in the figure. A consequence of using approximate 
vacua is that the equation of motion is violated over the entire domain of r^. 

From figure 11 we see that the symmetry breaking fluctuations grow and become non¬ 
linear on fixed s slices, exactly as expected for perturbations associated with the y direction. 
To check this is not a numerical artifact, we ran the same simulation using two different 
choices for the finite-difference stencil. In all cases the resulting evolution was very similar 
as seen in figure 13. This, coupled with the preservation of the symmetry for the exact in¬ 
stanton initial condition and the explanation of the linear instability in terms of parametric 
resonance [1], strongly indicates that the boost symmetry breaking is a real effect and not an 
artifact of the the discrete lattice spacing and time steps. Note that the strong breaking of 
the boost symmetry is absent in the solution to the dimensionally reduced dynamics (4.3). 
Thus, this is an effect which is only captured by intrinsically higher-dimensional simulations 
(in this case three-dimensional). 


5 Bubble Collisions with Fluctuations 

In section 4 we verified the accuracy of our lattice code and demonstrated the need for a 
full three-dimensional treatment of the two-bubble problem even in a simple case where we 
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Figure 10. The instanton profile for the linear symmetry breaking potential (2.1) with S = 0.1. We 
compare the exact instanton solution with various thin-wall approximations. On the left we show 
the numerically generated instanton profile used to initialize the bubbles. For comparison, we also 
include the thin-wall approximation, both in the case when we determine the vacua exactly (4.7), and 
in the case when we approximate the vacua to O(S^) (4.8). In the right panel we plot the violation of 
the instanton equation (2.3) for both our numerical profile and the thin-wall approximations. When 
measured by the maximal violation of the equation of motion, our solution is ten orders of magnitude 
more accurate than the thin-wall result. The mapping parameters were d = 0.35 and L = 1.6Ro. 
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Figure 11 . Field slice through the collision plane for a pair of colliding bubbles with initial conditions 
given by (4.8) in the linear symmetry breaking potential with 6 = 0.1. The rotational symmetry about 
the collision axis is maintained to high-precision, but one can clearly see the loss of boost invariance 
as the system evolves. The breaking of the boost symmetry develops on constant s slices exactly as 
one expects for resonantly amplified fluctuations with y dependence. 


break the boost symmetries but not the rotational symmetry. We now include a complete 
set of bulk fluctuations that break both the rotational symmetries and the boost symmetries. 
As we demonstrate shortly, the dynamical amplification of these fluctuations can lead to 
a complete breakdown of all of the boost and rotational symmetries in the system. This 
section constitutes the core part of our analysis. We consider several classes of single-field 
potentials, and we find that the extreme breaking of the spacetime symmetries is restricted 
to double-well potentials with mildly broken Z 2 symmetry. Section 6 extends the analysis to 
a two-field case. 

From our experience with the boost breaking fluctuations in section 4.3 and the linear 
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Figure 12. Left: The same view of the field as figure 8, except for an initial bubble profile (4.8) 
rather than the numerical result. Right: Contours of energy density — for the same simulation. 
A comparison with figure 8 shows that the collision results in the complete destruction of the boost 
symmetry. An animation of the field evolution is available at www.cita.utoronto.ca/~jbraden/ 
Movies/linear_delO.l_thinwall_field.avi. 



Figure 13. Direct comparison between our one-dimensional simulations using a 10th order Gauss- 
Legendre integrator and Fourier pseudospectral approximation for spatial derivatives, and our full 
three-dimensional parallel lattice simulations. We plot the value of the field at the initial point of 
collision between the two bubbles as a function of time. If we use the correct initial instanton profile 
the two results agree extremely well, apart from a small shift in the time of the initial collision. 
However, once we break the boost symmetry by using approximate choices for the locations of the 
true and false vacuum large deviations between the solutions appear. This effect is not captured by 
a one-dimensional lattice that assumes SO(2,l) symmetry for the evolving field. 
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analysis performed in [1], it is clear that asymmetric fluctuations will undergo resonant am¬ 
plification for certain types of bubble collisions. Eventually the amplified fluctuations will 
begin to interact nonlinearly. At this point the split between background and fluctuations 
becomes blurred and we must study the full three-dimensional nonlinear field theory. For a 
planar symmetric wall-antiwall pair, we demonstrated in [2] that the amplification of fluctu¬ 
ations transverse to the collision axis eventually lead to the dissolution of the walls and the 
creation of a population of oscillons distributed within a narrow transverse slab centered the 
site of the collision. Since we are considering bubbles here, the “transverse” wavenumbers 
can be split into a wavenumber associated with the radial hyperbolic direction y and one 
associated with rotation about the collision axis 9. For definitions of y and 9 see (4.1). The 
fluctuations with y dependence only were studied in section 4.3 where we looked at boost 
symmetry breaking fluctuations. In the remainder of the paper we explore the full evolution 
of three-dimensional bubble collisions with small initial fluctuations around the instanton 
profiles for a variety of potentials. In particular, we study the final outcome of the nonlinear 
interactions of the (linearly) amplified fluctuations. For cases where the bubble walls bounce 
off of each other many times, we will see that the results match our intuition developed from 
the planar wall limit studied in [2]. 

We take the initial bubble separation and the initial amplitude of the fluctuations to be 
independent free parameters. However, if we wish to study collisions between typical bubbles 
this will not be true. The action of the bounce, which determines the nucleation rate and 
therefore the typical bubble separation, scales as A“^. Meanwhile, the RMS amplitude of 
the fluctuations scales as V^. Therefore, increasing the amplitude of the fluctuations has 
the effect of increasing the nucleation rate and thus decreasing the typical bubble separation. 
Our primary motivation to treat these parameters independently was numerical. We simulate 
a finite volume with periodic boundary conditions, and therefore we have a finite amount of 
time before the bubbles begin to interact with their images and we are no longer simulating 
a simple two bubble collision. To study the effects of nonlinear interactions, the exponential 
growth of the fluctuations must push them into the nonlinear regime before the effects of 
periodicity reach the collision region. We chose the initial amplitude (while maintaining 
^ <C 1) to ensure this was the case. Of course, since the nucleation events are stochastic, 
for a given model a range of separations between nucleated bubbles will occur. 

5.1 Thin-Wall Double Well Case 

Our first example is the collision of two thin-walled bubbles in the linear symmetry breaking 
potential with 5 = 0.1. Aside from the new scale associated with the radius of the bubbles, 
this case is qualitatively the same as the collision of a pair of planar walls in the same 
potential. In section 4.3, we showed that perturbing the bounce solution slightly can lead to 
a dramatic breaking of the two boost symmetries. From our previous study of planar walls [2], 
we expect that the evolution with a full complement of initial fluctuations will ultimately 
lead to the dissolution of the bubble walls and the creation of a collection of oscillons in the 
collision region. Figure 14 demonstrates that this is indeed the case. Near the centre of the 
collision, the bouncing of the walls amplifies the (initially small) fluctuations. At mt ^ 40, 
these fluctuations become of similar amplitude to the oscillations of the SO(2,l) background. 
Shortly after this, the distinction between background and fluctuations in the collision region 
breaks down and the field rapidly condenses into a population of oscillons. In our choice of 
time coordinate, the condensation occurs first near the center of the collision. Meanwhile, the 
outward propagating torii produced during the first few collisions of the walls develop ripples 
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that eventually pinch off leading to the creation of “rings” of oscillons near the outer edges of 
the collision region. The two mechanisms described above correspond to the two production 
mechanisms we anticipated in section 4.1: parametric amplification of ffuctuations by the 
oscillating background field near the collision center, and the growth of ffuctuations on the 
outward propagating torii. 
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Figure 14. Development of the instability for two colliding bubbles in Minkowski space including bulk 
fluctuations corresponding to A = 10“^ around the thin-wall approximation. In the top row we plot 
the field distribution on a 2-D slice through the plane where the bubbles collide {horizontal projection) 
and on a slice parallel to the collision axis through the centers of the bubbles {vertical projection). In 
the bottom row we show contours of the energy density. Our lattice has N = 1024 points per side with 
a box size of mL = 13Ro = 13V26~^ and a time step dt = dx/12.3. The bubbles are nucleated with 
their centers separated by a distance mRsep = Videos corresponding to both sets of figures 

can be found at www.cita.utoronto.ca/~jbradeii/Movies/linear_delO.l_wfluc_field.avi and 
WWW. cita.utoronto.ca/~jbraden/Movies/linear_delO.l_rhocontours.mp4. 


5.2 Thin-Wall with Plateau 

For our next example potential, we modify the region into which the field tunnels by inserting 
a long-fiat plateau rather than a second well. This new potential is given by 

= 4>lf - + Cf>true 

\ Vq ~ e((/) — (ptrue) • 0 ^ 4^true 

where Vq = \{4>irue ~ ~ ^^o^true + Vq and (assuming 5 > 0) (ptrue is the largest solution 

to ~ ~ 0* This is meant as a toy example where the field tunnels out onto 

a fiat inflationary plateau, with the slope determined by 6 > 0. The portion of the potential 
traversed by the instanton as it tunnels is unchanged from the double well case (2.1), and 
therefore the initial bubble profile is the same. However, we expect the collision dynamics 
between this case and the thin-walled double-well to be radically different. For a sufficiently 
energetic collision, the free passage approximation [25, 62-64] will hold shortly after collision. 
This tells us that the collision will displace the field a distance Aj) ^ 200 down the plateau. 
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Unlike the double-well, this does not result in a restoring force pulling the field back towards 
the false vacuum. Therefore, there is no bouncing of the walls or other oscillations of the 
background field. Thus, there is no mechanism to pump fluctuations. Rather, the effect of 
the collision is to create a field gradient interpolating between the field value at the tunnel 
out location and its displaced value down the plateau. This gradient then propagates away 
from the collision as seen in figure 15. In front of the gradient, the field is at 0 (^true ~ 0o 
while behind it the field has been displaced down the plateau to (j) 

rue + 200 300. 

Although this behaviour is very similar to that of a propagating domain wall, the gradient 
does not interpolate between local minima of the potential. Thus, unlike a domain wall there 
are no bound state fluctuations associated with the interpolating field because V" is never 
negative. Somewhat suprisingly, the field at the collision site does not remain stationary 
on the plateau. Rather, it begins to slowly retreat back to the tunnel out location, at least 
for the duration of our simulations. From figure 17, we see that the field profiles both with 
and without initial fluctuations are nearly the same, indicating that the fluctuations play a 
subdominant role for this choice of potential. As a result, this is a situation in which the 
SO( 2 ,l) assumption provides an accurate description of the evolution. 

Now consider the effect of adding a nonzero slope to the flat plateau. Inside the bubble 
the field begins to roll down the plateau once it tunnels. Again, the collision effectively 
displaces the field further down the plateau and produces a steep gradient which propagates 
into the bulk. This gradient propagates on top of the previous rolling of the field down the 
potential. As with the flat plateau, figure 17 shows that the resulting evolution is insensitive 
to the choice of whether or not we include fluctuations. Therefore, we expect this situation 
to also be well described by 1 + 1 -dimensional simulations. 

5.3 Thick-Wall Double Well Case 

As our final potential we consider a collision between two thick-walled bubbles in the cubic 
symmetry breaking potential with 6 = 0.99. In figure 18 we show the result of one such 
collision, where we have made the assumption of SO( 2 ,l) symmetry in order to run a very 
high resolution one-dimensional simulation. The collision still leads to large oscillations of 
the field around the minimum, but the the field no longer becomes temporarily trapped in 
the false vacuum minimum between collisions. Instead of having repeated collisions between 
a pair of walls, we instead have oscillations of the field around the true vacuum. Within a 
few oscillations, the amplitude damps to < O.50o- We already know from our study of single 
field bubbles that spatially homogeneous oscillations of the field with this amplitude do not 
lead to a strong instability. Since the oscillations from the collision have additional spatial 
localization, it is clear that the fluctuations will not experience significant amplification in 
this case either. We see this directly in figure 19, where the rotational symmetry about the 
collision axis is preserved to a high degree during the collision. As well, figure 20 further 
shows that the boost symmetry is extremely well preserved. 

5.4 Bubble-Domain Wall Collision 

Thus far we have only considered the collision of two bubbles in a reference frame in which 
they nucleate at the same time. We will now address the opposite regime where one of the 
bubbles is much larger than the other. To model this situation we consider the collision of a 
bubble with a planar domain wall, with the planar wall meant to be a substitute for the large 
bubble. In an actual collision between false vacuum bubbles, the wall has a very large Lorentz 
contraction factor 7 1. This creates a very large hierarchy between the size of the smaller 
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Figure 15. A collision of two single-field bubbles with a flat plateau onto which the field tunnels (5.2). 
The horizontal field slices are orthogonal to the collision axis and through the center of the collision 
region. For orientation, the back panel is a slice parallel to the collision axis, illustrating the growth 
of the two bubbles. Blue corresponds to regions where the field is near the false vacuum, red regions 
where it is near the tunnel out point, and green regions down the plateau away from the tunnel out 
point. A illustrative video of the process can be found at www. cita.utoronto. ca/~ jbraden/Movies/ 
collision_plateau_delO.l_field.avi. 


bubble and the thickness of the planar wall which we are unable to resolve with a fixed lattice. 
To avoid this problem we restrict ourselves to mildly relativistic walls. A physical situation 
where this occurs is finite-temperature phase transitions, where interactions between the 
bubble wall and the surrounding thermal medium prevent the speed of the bubble walls from 
nearing the speed of light. Of course the initial spectrum of fluctuations is different in the 
thermal case and interactions with the surrounding thermal bath must also be accounted for. 
Therefore a direct comparison with our results cannot be made, although the same qualitative 
behaviour will persist. As well, we consider only collisions in the linear symmetry breaking 
potential (2.1), again with 6 = 0.1. From our results for the collision of two bubbles nucleated 
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Figure 16. A collision of two single-field bubbles with a tilted plateau e = 0.01. The field slices are 
the same as in figure 15. Inside the bubble, the field now rolls down the potential. However, the result 
of the collision is essentially the same as in the case e = 0. The field is displaced down the potential 
at the collision, and a field gradient then propagates into the interior of the bubble. Aside from this, 
the field rolls down the potential due to the constant force resulting from the slope. 


at rest at the same time, it is clear that these types of collisions will again result in strong 
amplification of the fluctuations and the eventual production of a population of oscillons. 
The only new ingredient here is that the planar wall and bubble wall no longer carry an 
equal and opposite amount of field momentum density along the collision axis = (j)d^\\(j). 
Therefore, if we consider a small box around any region where the planar wall and bubble 
are undergoing a collision, the energy within that box carries a nonzero momentum along the 
collision axis. As well, the collision no longer occurs in a single spatial plane, but rather on 
a curved hypersurface defined by all of the instantaneous intersections between the bubble 
and the planar wall. As a result, the collision products are produced along a nonplanar 
surface with a nonzero velocity (both relative to our simulation coordinates). Since we are 
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Figure 17. Comparison of bubble profiles orthogonal to the collision axis for a simulation with initial 
fluctuations Sc/) ^ 0 {black line) and without initial fluctuations Sc/) = 0 {red triangles). The choice 
of potential and initial conditions are the same as figure 15 {left) and figure 16 {right). We plot four 
times mt = 34.47, mt = 45.96, mt = 57.45 and mt = 80.43. The width of the region displace from 
the false vacuum (0 = — ^q) increases with time. The two profiles match extremely well at all times, 
indicating that the fluctuations are stable for this class of bubble collisions. The oscillations near the 
edge of collision region are numerical artifacts. In the left panel we show the result for a fiat plateau 
e = 0, while on the right we have e = 0.01. For the purposes of plotting, we have downsampled our 
output grid by a factor of 2. 
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Figure 18. The evolution of a pair of colliding bubbles in the cubic potential with 5 = 0.99 under 
the assumption of exact SO(2,l) symmetry. We plot the field as a function of time and position along 
the collision axis. 


not properly boosting an initial two-bubble field configuration to set the initial conditions, 
the collision hypersurface will be nonplanar in any reference frame. 

Figure 21 shows that these expectations are true. As the wall accelerates into the false 
vacuum it partially engulfs the expanding bubble. As with the two bubble case, the field 
experiences a series of excursions back to the false vacuum. These drive an instability in 
linear (non rotationally invariant) fluctuations analogous to the instability in the pair of 
bubbles. The torii resulting from these excursions no longer distribute themselves in a single 
plane, but rather they move in the same direction as the wall. As a result, the oscillons that 
are ultimately produced from these growing fluctuations are distributed in some thin curved 
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Figure 19. A collision of two-thick wall bubbles in the potential (2.2) with S = 0.99. The field slices 
and coloring are the same as in figure 14. A video of the field evolution is available at www.cita. 
utoronto.ca/~jbraden/Movies/thick_wall_collision.avi. 


hypersurface rather than a plane. 

6 Extension to Two Fields 

By considering the simple case of a single-field potential, we discovered that the dynamics of 
bubble collisions can be considerably more intricate than expectations based on 50(2,1) or 
50(2) symmetry for the field profile. In particular, for single-field double-well potentials with 
moderately broken Z 2 symmetry we found that the collision of a pair of bubbles ultimately 
leads to the production of a population of oscillons. These types of symmetry breaking 
potentials are often used to model first-order phase transitions, and thus our results may 
have some relevance in that context. However, the bubbles must coalesce in order for the 
phase transition to complete. When the bubbles begin to coalesce we must consider collisions 
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Figure 20. Evolution on a line orthogonal to collision axis through the initial point of collision for a 
collision between thick-walled bubbles in the cubic symmetry breaking potential with S = 0.99. We 
have taken a fairly large initial fluctuation amplitude to demonstrate the preservation of the boost 
symmetry. The fluctuations were those in the Minkowski false vacuum with A = 10“^. 


between many bubbles rather than just two, leading to an explicit and large breaking of the 
SO(2,l) symmetry. 

A natural scenario to consider collisions between isolated pairs of bubbles is in false 
vacuum eternal inflation. Unfortunately, it seems exceedingly difficult to both realize inflation 
and produce oscillons in bubble collisions using a single-field potential with only two minima. 
The simplest way to have a viable single field open inflationary model produced by the 
nucleation of a bubble is to have the field tunnel out onto a plateau that can support A In a ^ 
50 efolds of inflation.^ Our mechanism for oscillon production relies on oscillations of the field 
around the true vacuum within the collision region. However, if the plateau is long enough 
to support a sufficient period of inflation, any rolling motion of the field along the plateau 
must experience many efolds worth of Hubble overdamping. Therefore, the requirement of an 
inflationary plateau makes it difficult for the field to rebound off of a “wall” in the potential 
and return to the false vacuum side of the barrier shortly after collision. Thus the basic 
mechanism by which initial fluctuations are amplified disappears. 

However, embedding inflationary models based on bubble nucleation into a realistic 
high-energy theory will likely involve considering models with many scalar fields.^ With this 
in mind, consider a simple two-field potential 


V{a,4>) = ^{a^-aif + X^S 


a^ao 




- + — + Act— (<7 - aoY 


+ Ao-CqC*/) + Vq . 

(6.1) 


Schematically, the potential has the form VtunneM) + Vcoupling{<^,(t>) + Vinfiation{4>)^ where 
we have chosen the particular form Vinfiationi(p) = KcrQecf) + Vo- Our choice Vinfiation = 
\cr(yQe(f) + Vb can be viewed as a linearization of the potential around the tunnel out region. 
There is a local minimum at cr —(Tq and 0 « 0, and a long trough at <7 uo along which a 
is heavy and (j) is light. Since we now have two mass scales at our disposal, we can accomodate 


^Another possibility is that collision displaces the field to a location where it begins to slow-roll, thus 
starting inflation. 

^Even this statement is probably too simplistic, as the correct high-energy theory may not even be de- 
scribable in terms of a scalar field theory. 
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Figure 21. Field and energy density evolution for a collision between a planar domain wall and a 
single nucleated bubble. In the top row we show two slices of the field configuration. The first slice 
(projected vertically) is perpendicular to the planar wall and cuts through the center of the bubble. 
The second slice (projected horizontally) is parallel to the domain wall and displaced slightly from 
the center of the bubble towards the initial location of the domain wall. As the wall sweeps past the 
bubble, the instantaneous collision location moves along the collision axis, as do the tubes of false 
vacuum and oscillons produced by the collision. In the bottom row, we instead show contours of 
energy density for the same three time slices. 


the production of oscillons by exciting the a field while simultaneously permitting slow-roll 
inflation along the (j) direction. By adjusting Vinfiation we can effectively reproduce any model 
of single-field inflation we wish. 

We only consider parameter choices such that the tunnelling dynamics is dominated 
by a, while the subsequent post-tunnelling evolution is dominated by (/). Although we have 
not performed an exhaustive study for different choices of when = 1 this behaviour is 
generic for the thin-walled case d <C 1. This is illustrated in figure 23, where we explicitly 
see that during tunnelling the field first moves almost exclusively in the a direction before 
making a sharp turn at the end so that it is moving along the slow-roll plateau. Effectively, 
the behaviour of the tunnelling field a is fixed by the double well structure, and the inflaton 
field (j) simply reacts to the presence of the domain wall in a. In figure 24 we show the 
potential seen by the effective field G?Xe// ~ tunnels through the barrier, as 

well as a comparison to an analogous single-field potential. The structure of the potential as 
seen by the field while it tunnels and subsequently slow-rolls is very similar to the double-well 
with an appended slow-roll plateau (5.2). As well, notice that there is no steepening of the 
potential near the beginning of inflation as is often assumed in single-field models of open 
inflation from bubble nucleation. 
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Figure 22. The two-field potential (6.1) for e = 0.01, S = 0.2, = 1 and Vq = 0. We have in 

mind a scenario where the field is initially trapped near the local minimum at a ~ —cro,0 ~ 0 and 
subsequently tunnels into the nearly flat trough at a = cro- For clarity we have clipped the potential 
for V > A(7^. 



Figure 23. The instanton solution used to set initial conditions for our numerical simulation. In 
the left panel we show the path of the instanton in field space superimposed on isocontours of the 

potential. In the middle panel we show the profile of the tunnelling field a and the “inflaton” field 0 

as a function of the Euclidean radius r^;. In the right panel we plot the spectral coefficients for each 
of the fields, from which we can see the exponential convergence of the series and the roundoff plateau 
arising around i ^ 80. The model parameters were (5 = 1/5, e = 0.01 and = 1. Our mapping 

parameters were L = and d = 0.5, and we used 100 mode functions. 


Figure 25 and figure 26 illustrate the collision dynamics between two of these multifield 
bubbles. We do not perform an exhaustive study of possible behaviours as a function of 
initial bubble separation or model parameters. Instead, we simply choose a set of parameters 
and initial conditions to demonstrate that oscillons can be produced in bubble collisions in 
this model. For our choice of parameters, the collision causes a large excitation in a and its 
subsequent evolution is very similar to the single-field double well case. However, for a fixed 
a cur 7 ^ ^0 the potential V {(Tcur^ 0) has a minimum at 0 = —ea^/{a — (Jq)^ <C o-q. As a result, 
the oscillations of a pull the field (j) back towards the origin, thus undoing the previous 
rolling down the plateau. The evolution of the a field amplifies fluctuations and leads to 
the creation of oscillons in the a direction. This is analogous to the single-field evolution. 
However, in the cores of the oscillons, the a field makes large excursions away from (Jq. As 
a result, within these cores (j) remains trapped at the origin while the field outside of these 
cores once continues to roll down the plateau. In the inflationary setting, there will still be a 
vacuum energy Vq sufficient to drive a period of slow-roll inflation, so as the oscillons dilute 
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Figure 24. The potential seen by the field as it tunnels in the Euclidean radial direction {blue line). 
Here we have defined dxlff = da^ + dcj)^ as the path length in field space. For comparison, we also 
include the potential E(0, Xeff) that would be seen by the field if we had instead tunnelled a distance 
Xeff with (j) = 0 {red dots). In this latter case, the potential is then the same as our single-field cubic 
symmetry breaking potential. 
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Figure 25. Evolution of the “infiaton” field (p in our two-field model with e = 0.01 and 6 = 0.2. The 
horizontal projection is a slice along the collision axis through the center of the bubbles and the vertical 
slice is orthogonal to the collision and centered on the collision site. Red corresponds to 0 values near 
the origin (i.e. in the false vacuum or at the tunnel out location). The pips where the field is pulled back 
up the potential are locations where the a field is fracturing into oscillons. An animation corresponding 
to this evolution is available at www. cita.utoronto. ca/~jbraden/Movies/twof ield_inf laton. avi. 


we expect inflation to eventually restart. A novel feature of this setup is that (j) will be quite 
inhomogeneous (within the collision region) at the start of this new inflationary phase, and 
we expect this inhomogeneity to persist until the oscillons decay. 

6.1 Inflationary Model Building and Consistency of the Minkowski Approxima¬ 
tion 

Let’s now consider the consistency of our approach and application of this type of model to 
inflation. In our numerics, we have assumed that the instanton profile is well approximated 
by taking the background to be Minkowki. As well, for our lattice simulations we have 
further assumed that we can account for expansion of the universe by taking a uniform 
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Figure 26. Evolution of the “tunnelling” field a in our two-field model with e = 0.01 and d = 0.2. 
As with previous figures, blue corresponds to a ^ —cfq (i.e. near the false vacuum) and red to a ^ cfq 
(i.e. in the “inflationary” trough). The vertical slice is parallel to the collision axis and the horizontal 
slice is orthogonal to the collision axis. An animation corresponding to this evolution is available at 
WWW. cita.utoronto.ca/~jbraden/Movies/twofield_tunnel.avi. 

fixed Hubble constant.Relaxing this latter assumption is a significant numerical challenge 
and we therefore leave it to future work. This is required to obtain accurate predictions 
for the observational signatures that remain after inflation has occurred within the bubble. 
While this is certainly an interesting question, in this paper we are primarily concerned with 
the dynamics of the collision itself. The relevant time scale is then (ActCTq) ^ , and we 
assume this is much shorter than a Hubble time in the ambient spacetime in order for our 
approximations to be valid. 

First consider the restrictions imposed by our computation of the instanton profiles. 
The least stringent requirement is that the CdL instanton exists. Roughly, this requires that 
the bubble wall fit within a Hubble radius determined at the local maximum of the potential 
^max' For our setup, the thickness of the wall is determined by Vaaio- = 0 , 0 = 0 ) ^ ActCTq, 
and therefore we require ActCTq ^ Vo/3Mp + AcrCrQ/12Mp =4> 12^^ ^ 1 + A 

more stringent constraint comes from requiring that the initial radius of the bubble is much 
less than the Hubble scale in the ambient spacetime We have rrieffRinit ^ and 

HjyMp ^ + ^, so this gives ^1 + 4 ^^) < S. In the limit that Vq > 5A^(Jq, 

this gives ^ ^ 2^-2 For the opposite limit we instead get <C S. 

Finally, consider the restriction imposed by assuming the approximation of a fixed 
homogeneous background Hubble. This constraint only needs to be fulfilled if we want an 
approximate description of the dynamics prior to the collision, or if we want to track the 
long time evolution after the collision. The average expansion rate inside and outside the 
bubble must be much greater than the difference in these two expansion rates. As well, for 
the case when the field slow-rolls inside of the bubble, our approximation will only be valid 
when the difference between the Hubble parameter in the center of the bubble (where the 

^°The figures presented in this paper were taken from simulations using a Minkowski background. However, 
we also did several runs with a fixed Hubble constant H over the entire simulation volume. The main effect 
of this expansion was to delay the time that the bubbles first collide, while the collision dynamics itself was 
unaffected. 
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field has rolled furthest) and the edge is again much less than the average. This then requires 
Vq ^ ^ActCTq, |AcrCrQe0x=o(^)| where (/)x=o is the value of 0 at the center of the bubble at time 

t. 

7 Observational Prospects 

In the previous sections we demonstrated the presence of a previously ignored instability 
in the fluctuations around colliding vacuum bubbles in double-well potentials. As a result 
of these instabilities, the near SO(2,l) symmetry of the initial configuration becomes badly 
broken during the course of the dynamical evolution. This eventually results in the production 
of a population of oscillons in the collision region. In this section we briefly comment on some 
possible implications of these results, restricting ourselves to potential signals that rely of 
the breakdown of the SO(2,l) symmetry. This means we will not discuss possibilities such 
as a large scale modulation of coupling constants. Modulations could result from single- 
field collisions in potentials with a plateau if the vev of the field forming the bubbles sets 
an effective coupling constant in the full theory. There are many scenarios in which one 
can imagine embedding our mechanism, and we will distinguish between two cases that we 
refer to as super horizon and subhorizon collisions. In the first, we assume that our entire 
observable universe fits within one of the nucleated bubbles, with the nucleated bubble itself 
embedded in some parent false vacuum. The collision is then with a neighbouring bubble 
which has also nucleated within the parent vacuum. Of course, this is the standard scenario 
for testing open inflation following bubble nucleation in an ambient false vacuum. Past 
studies of signatures from such collisions are based on the assumption of SO(2,l) symmetry, 
so in this case in particular we would like to address any novel implications of the breakdown 
of the symmetry. One could also imagine scenarios where our observable universe forms in 
the future light cone of the collision, although we will not discuss this possibility here. In the 
second, we instead assume that the bubble nucleations are occurring on subhorizon scales 
within our Hubble volume. These nucleations could either have occurred in the past (such 
as an early first order phase transition within our Hubble volume) or in the present. The 
reader should keep in mind that we have not performed a detailed analysis for any of the 
speculated possibilities in this section, so the magnitude of many of these effects may prove 
to be undetectably small. We plan to provide a more detailed study in a future publication, 
including the effects of gravitation as required to obtain the present day signal. 

7.1 Superhorizon Bubble Collisions 

First consider the case where our entire observable universe is contained within a nucleated 
bubble. Our results demonstrate rather clearly that observational signatures of collisions 
do not necessarily posess an azimuthal symmetry about the collision axis as is assumed in 
the existing literature. A remnant of the azimuthal symmetry persists in the sense that the 
effects of the collision are confined to a circular disk of the sky, but the interior of this disk 
need not posess any additional symmetries. 

Perhaps the most interesting possibility is the production of gravitational waves by the 
fracturing of the bubble walls at the onset of nonlinearity amongst the fluctuations. This 
effect is absent in an exact SO(2,1) collision, since the symmetries forbid the presence of tensor 
modes. Initially the gravitational waves are produced with subhorizon sized wavelengths, 
but as the oscillons dilute and inflation restarts they are stretched outside the horizon. The 
resulting signal will have a characteristic wavelength as well as a directional dependence on 
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the sky that may be detectable in polarization data. As well, the directional dependence of 
this polarization signal will be highly correlated with other possible signatures, for example 
a hot or cold spot produced by the collision. Since the amplitude of the waves will damp 
until they are stretched outside of the horizon, the magnitude of the effect will depend on 
the typical scale of the nonlinear fracturing of the wall relative to the Hubble scale inside the 
bubble. 

Another effect is related to sign of the temperature perturbation induced by the collision. 
In the single-field case with the field tunnelling onto a plateau, the collision simply displaces 
the field down the plateau. Therefore there are fewer efolds of infiation in the collision region 
than the rest of the bubble, so 51n(a)|p < 0 and we obtain a hot spot. This was verified 
in [28] for a specific choice of infiat ionary model. However, when considering more complex 
models this simple intuition will not hold in general. For example, in our two-field infiationary 
model (6.1) the “infiaton” (f) within each bubble has already begun to roll down the trough 
prior to the collision. The effect of the collision, at least for the parameters used in this paper, 
is to pull (j) back up the potential and to excite oscillations in a. Eventually the dynamical 
evolution of a produces oscillons, with the novel feature that the field (j) is trapped Sit (/) ^ 0 
in the interior of each oscillon. If the net effect of the collision were to simply pull (j) back 
up the potential, this would prolong the infiationary epoch and produce a positive 5\n{a)\p. 
However, the full dynamics creates additional contributions to the energy density (initially 
the oscillations in a and later the oscillons). Thus, a detailed simulation including the effects 
of gravity is probably needed to ascertain the final sign of 51n(a), although it seems likely 
that positive values for the perturbation and thus cold spots can be generated. 

Finally, in our two-field model the infiationary epoch within the collision region restarts 
from a highly inhomogeneous state due to the population of oscillons. As long as these 
oscillons persist, they trap the infiaton field at the origin in regions of size = (ActCTq)”^/^. 
Whether or not any remnant of this (non-vacuum) initial state persists is a very interesting 
question. One possibility is that that once the oscillons decay, the pips of infiaton field 
that were held near the origin will be stretched in physical size as the universe expands and 
could end up sourcing curvature fiuctuations As well, the field fiuctuations in the collision 
regions will have initial conditions for the subsequent infiationary era that non-vacuum and 
highly nonGaussian. If some trace of these infiaton perturbations survives, it would provide 
an example of spatially localized nonGaussian density fiuctuations. 

All of the effects speculated above will be rather strongly constrained by data, so a cer¬ 
tain amount of tuning must be applied in order to construct models that produce signatures 
that are possibly observable but not yet ruled out by data. As well, additional tuning is 
needed in order for the underlying theory to predict a nonnegligible number of potentially 
observable collisions for a typical observer. Our purpose in this paper was simply to un¬ 
derstand the dynamics of individual collisions so we will not touch on these model building 
issues here. However, they must be addressed in order to assess the probability of seeing a 
given signal and to ultimately constrain the underlying theoretical model. 

7.2 Subhorizon Bubble Collisions 

Now let’s consider the case when the bubble collisions occur in some first-order phase transi¬ 
tion within our horizon. Such a transition may have occurred in the past when the universe 
was extremely hot. However, we can also envision a scenario where the phase transition is 
happening around the present time. An interesting possibility would be for the field respon¬ 
sible for dark energy to undergo a first order phase transition. As in the previous subsection. 
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we restrict ourselves to signatures that result from collisions and the breaking of SO(2,1) sym¬ 
metry, although there are additional possibilities than are not of this type. The discussion 
here in many ways mirrors the discussion for the superhorizon collisions, with the possible 
signals being more spatially homogeneous counterparts to the highly spatially localized effects 
discussed in section 7.1. 

First consider the production of gravitational waves during a first order phase transition. 
When the bubbles collide, gravitational waves are produced with wavelengths of order the 
typical size of the bubbles at collision. This source of gravitational radiation is well known. 
It is often studied using the envelope approximation, which neglects all of the field collision 
dynamics and treats the bubbles as if they were undergoing free expansion. However, if the 
phase transition is not too rapid (so that the bubbles have a chance to grow before colliding), 
then we expect that the fracturing process explored in this paper will occur in individual 
collisions. This will lead to an additional peak in the gravitational wave spectrum associated 
with the size of the unstable modes produced during the collision. These two sources of 
gravitational waves mentioned above are a direct result of scalar field dynamics and are 
present even in vacuum. However, in a high temperature phase transition the subsequent 
(turbulent) dynamics of the plasma may source additional gravitational waves. In this case, 
the fracturing of the walls and long-lived oscillons may force turbulent motion and thus 
indirectly produce additional gravitational waves. 

The production of oscillons during the phase transition also perturbs the expansion 
history of the universe by adding a component that dilutes as collisionless dust to the ambient 
cosmological constant (vacuum transition) or radiation bath (high temperature transition). 
If the oscillons form a significant fraction of the post-transition energy density, this could 
lead to a temporary stage of matter domination in the early universe. Since the linear 
growth of subhorizon adiabatic density perturbations is controlled by the expansion history, 
a temporary stage of matter domination could leave an observable remnant in large scale 
structure today. 

8 Applications to Other Scenarios 

We have restricted ourselves in this paper to some very simple toy potentials that permit 
Coleman-deLuccia bubble nucleation. However, the basic mechanism — amplification of non 
SO(2,1) symmetric fluctuations — likely plays a key role in many other scenarios. For a 
small number of fields, the important feature of the potential appears to be the presence of 
a slightly asymmetric double-well structure in one of the field directions. In this section we 
provide several examples from the recent literature where the asymmetric fluctuations were 
(we believe incorrectly) neglected. 

A direct application is to the bubble baryogenesis scenario of [65]. This scenario takes 
a two-field model with slightly broken internal global 0(2) symmetry. Baryon number is 
generated by the nucleation and expansion of bubbles, and also during bubble collisions. 
The collisions considered by the authors of [65] very closely resemble our results for the thin- 
walled bubbles, except that S0(2,l) symmetry was explicitly enforced. It would be interesting 
to see how the net baryon number production is influenced by breaking of S0(2,l), as well 
as the extent to which baryon number is trapped within the oscillons. 

Another possible application of these results is to boom and bust inflation [66] and 
other models of flux discharge cascades driven by bubble nucleation and expansion in a 
compact extra dimension [67, 68]. However, there are several caveats to a direct application 
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of our results. First of all, the bubble walls in these models can move at extremely relativistic 
speeds and are thus extremely Lorentz contracted in the center of mass frame for the collision. 
Therefore, the duration of each collision is very short and the walls may simply pass through 
each other before the fracturing process has time to occur. A second caveat is that the 
bubble walls in the flux cascade picture are D-branes instead of scalar field kinks. Since the 
nonlinear interactions that lead to the breakup of the walls are sensitive to the details of the 
high energy completion of the theory, the nonlinear fracturing stage may be altered in the 
flux discharge case. 

A final recent scenario appearing in the literature where non-symmetric fluctuations are 
likely to play a key role but were not considered is in [69]. Here, the authors study late-time 
violation of the free-passage approximation by displacing the field to very near local maxi¬ 
mum of the potential in the collision region. However, they make the assumption of planar 
symmetry and work in 1 -h 1 dimensions. The nonplanar fluctuation modes will experience 
a tachyonic instability [60, 70] that quickly invalidates the planar symmetry approximation. 
In fact, the seeds of this are already in the author’s work as their effect relies on precisely 
this tachyonic instability for the k± = 0 mode. 

Beyond these obvious applications of the ideas presented in this paper, there are many 
other scenarios in which corrections to SO(2,l) symmetry may be expected. In this paper, we 
only considered collisions between bubbles formed from the same instanton. Therefore, there 
were no topological constraints preventing the walls of the colliding bubbles from eventually 
annihilating each other. A natural extension is to consider more complex potentials with 
many local minima. Multiple instantons describing different decay channels for a single 
metastable minima can then exist, and collisions can occur between bubbles with different 
interiors. In this case, the field values in the bubble interiors enforce topological constraints 
on the collision dynamics. Therefore, in general new domain walls interpolating between the 
bubble interiors will form in the collision region. These newly formed domain walls can have 
their own internal dynamics which can again drive instabilities that lead to a breakdown of 
the SO(2,1) invariance. This is the generalization of the oscillations of the planar shape mode 
in the double well potential studied in [2] to the case of bubble collisions. Since the internal 
dynamics generally emit radiation, and this radiation will not in general respect the SO(2,l) 
symmetry, any signatures that may result from it again differ from the results of an SO(2,l) 
simulation. 

9 Conclusion 

We performed fully three-dimensional simulations of bubble collisions in scalar field theories. 
Our treatment is novel because we include, for the first time, the effects of quantum fluctu¬ 
ations (in the semiclassical wave limit) on the dynamics. Recent interest in this topic has 
been driven largely by false vacuum eternal inflation and early universe phase transitions. In 
these instances, bubbles initially form through either quantum or thermal nucleation then 
subsequently undergo collisions. In the case of bubbles nucleating in a false vacuum, previous 
studies have assumed that collisions between two such bubbles possesses an exact SO(2,l) 
symmetry. This symmetry derives from the partial breaking of the SO(4) symmetry of the 
instanton solution. However, it is important to keep in mind that the instanton only describes 
the most likely bubble to nucleate. In reality, the actual shape of the nucleated bubble will be 
slightly deformed from the perfect instanton profile due to the quantum nature of the nucle¬ 
ation. Even ignoring these deformations to the initial bubble shape, the bulk fluctuations in 
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the ambient spacetime and inside the bubble are present after the nucleation of a bubble. In 
fact, these sub horizon fluctuations inside the bubble must be present for inflation inside the 
bubble to seed density perturbations for the standard hot big bang. In other applications, 
such as high temperature phase transitions, the nucleation of individual bubbles is again 
driven by coherent structures arising from stochastic fluctuations. Hence in all these cases 
the exclusion of fluctuations, which is implicit in imposing exact SO(2,l) invariance, is not 
consistent with the process of nucleation or creating a viable inflationary model. Since the 
fluctuations do not obey the assumed SO(2,l) symmetry (or SO(2) symmetry in the case of 
thermal nucleation), it is important to test these assumptions.^^ 

We studied collisions between pairs of nucleated bubbles in a variety of single-field and 
two-field potentials using highly accurate nonlinear three-dimensional lattice simulations. A 
novel aspect of this investigation was the development of a pseudospectral approach for find¬ 
ing SO(4) symmetric instanton solutions, rather than the ubiquitous overshoot-undershoot 
method. The accuracy of our instanton profiles are only limited by machine precision round¬ 
off errors, and the procedure easily generalizes to multifield potentials. Although we only 
considered the case of instantons in Minkowski space, our method can be trivially generalized 
to the case of a fixed background geometry, and with a little extra effort to the case of a 
dynamical metric coupled to the scalar fields. 

We studied two types of single-field potentials with radically different collision dynamics: 
double-well potentials with a broken Z 2 symmetry, and potentials with a single local mini¬ 
mum and a linear plateau which the field tunnels onto. Under the assumption of SO(2,l) 
symmetry, collisions in the double-well potentials with mildly broken Z 2 symmetry cause the 
bubble walls to bounce off of each other many times during the collision. These bounces 
produce an outgoing pattern of toroidal waves centered on location of the initial collision and 
expanding in the plane orthogonal to the collision axis. When we include symmetry breaking 
fluctuations to model quantum effects, the fluctuations experience a strong instability as a 
result of the bouncing motion in the SO(2,l) background solution. These fluctuations are 
quickly driven into the nonlinear regime, at which point a split into an SO(2,l) symmetric 
background and nonsymmetric fluctuations loses its utility. A fully three-dimensional non¬ 
linear description is required to describe the subsequent evolution. The nonlinear dynamics 
causes the bubble walls and expanding torii to fracture into a network of localized blobs 
of field known as oscillons. At this point, the SO(2,l) symmetry is spoiled completely and 
any dimensional reduction based on it will be a very poor approximation to the collision dy¬ 
namics. To the best of our knowledge, this is the first example where the three-dimensional 
nature of the problem plays an important role in the field dynamics for the collision between 
a pair of bubbles.Meanwhile, in situations where the field tunnels out far away from a 
minimum of the potential (either because the Z 2 symmetry of the double-well is strongly 
broken or there is a long plateau), we found no evidence that the SO(2,l) breaking fluctua¬ 
tions experienced instabilities.^^ As a result, these situations can be well-approximated by 
dimensionally reduced simulations. 

In the single field models we considered, it is difficult to have both an inflationary stage 

noted before, by this we mean that individual realization of the fluctuations do not preserve the 
symmetry, even if they do preserve the symmetry in a statistical sense. 

^^Easther et. al also performed three-dimensional bubble simulations [25], but they did not include fluc¬ 
tuations in their initial conditions and the effect that they discuss can be captured using 1-h 1-dimensional 
simulations. 

^^For the badly broken Z 2 double-well, this statement may be an artifact of our choice of potential, but for 
the plateau in the tunnel out region it should be generic. 
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in the interior of the bubble while simultaneously having collisions that drive an instability 
in the symmetry breaking fluctuations. This essentially follows from the fact that inflation 
requires a long region of the potential whose characteristic mass scale is much less than 
the Hubble constant. In contrast, the bouncing of the walls and subsequent production of 
oscillons only occurs if the typical mass scale of the potential minimum is much greater 
than the Hubble constant. For the case of identical bubbles, which was the focus of this 
paper, if the field repeatedly returns to the false vacuum in the collision region it must also 
repeatedly cross the portion of the potential that drives inflation. Since the contraints of 
the curvature of the potential for slow-roll inflation and oscillon production are mutually 
exclusive, this suggests it is extremely difficult to both drive inflation within the bubble 
and produce oscillons in collisions for typical single field potentials. To avoid this obstacle, 
we considered a very simple two-field model where we allowed one of the field directions to 
dominate the tunnelling and the other direction to drive slow-roll inflation. Since current 
theoretical ideas suggest there are many effective degrees of freedom at the energy scales 
relevant for inflation, the introduction of an additional field is quite natural from a top down 
perspective. During a collision between two bubbles in the two-field model, the evolution in 
the tunnelling direction closely resembles that of the single-field collision in the double-well 
with mild breaking. As a result, oscillons are able to form in this direction. The coupling 
between the tunnelling field and the inflaton direction then causes the “inflaton” field to be 
pulled back up the potential. This entire process again badly breaks the SO(2,l) symmetry 
and requires more than a one-dimensional simulation to capture. As well, the overall effect 
on the expansion history in the collision region is not necessarily to reduce the duration of 
inflation. This is in contrast to the single field case with the field rolling down a plateau, 
where the collision simply displaces the field down the plateau. Even in this relatively simple 
extension beyond single field models, the full collision dynamics is not well approximated by 
the theory of a single effective field. Therefore some caution should be used when drawing 
conclusions about collisions in models with a large number of fields. Since the fields now 
have many more directions to explore, both the tunnelling process and collision dynamics 
could change substantially in this limit. 

Our simulations are not fully relativistic, so we are unable to track the subsequent 
evolution of the fields through the complete phase of inflation inside the bubble. This is 
certainly an interesting question, but it is also a very challenging numerical problem that is 
much more difficult than the corresponding symmetry reduced (l-hl)-dimensional problem. 
However, since we can construct models where the time scale for oscillon formation is much 
less than the ambient Hubble time, our results set the initial conditions for the subsequent 
inflationary epoch within the collision region. Since the SO(2,l) symmetry is broken so 
strongly in some cases, this question cannot be addressed within the framework of a symmetry 
reduced 1-h 1-dimensional problem. In particular, previous signatures associated with bubble 
collisions have assumed an azimuthal symmetry will hold, which results from the assumption 
of SO(2,l) dynamics for the spacetime. 
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